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Abstract. Columnar defects provide effective pinning centers for magnetic flux 
lines in high-T c superconductors. Utilizing a mapping of the statistical mechanics 
of directed lines to the quantum mechanics of two-dimensional bosons, one expects 
an entangled flux liquid phase at high temperatures, separated by a second-order 
localization transition from a low-temperature "Bose glass" phase with infinite 
tilt modulus. Recent decoration experiments have demonstrated that below the 
matching field the repulsive forces between the vortices may be sufficiently large to 
produce strong spatial correlations in the Bose glass. This is confirmed by numerical 
simulations, and a remarkably wide soft "Coulomb gap" at the chemical potential is 
found in the distribution of pinning energies. At low currents, the dominant trans- 
port mechanism in the Bose glass phase proceeds via the formation of double kinks 
between not necessarily adjacent columnar pins, similar to variable-range hopping 
in disordered semiconductors. The strong correlation effects originating in the long- 
range vortex interactions drastically reduce variable-range hopping transport. 



1 Pinning of Flux Lines to Columnar Defects 

The remarkably rich phase diagram of magnetic flux lines in high-T c super- 
conductors, especially when subject to point and / or extended disorder, has 
attracted considerable experimental and theoretical interest (for a recent re- 
view, see Ref. [1]). For the purpose of applying type-II superconductors in 
external magnetic fields, an effective vortex pinning mechanism is essential, 
in order to minimize dissipative losses caused by the Lorentz-force induced 
movement of flux lines across the sample. This issue becomes even more 
important for the high-T c superconductors, because the strongly enhanced 
thermal fluctuations in these materials render the Abrikosov flux lattice un- 
stable, and one therefore expects a first-order melting transition leading to a 
normal-conducting flux liquid phase with well below the upper critical field 
H C2 (T), at least for ideally pure systems [2]. It is very difficult to pin such 
a vortex liquid consisting of fluctuating flexible lines by just intrinsic point 
disorder (e.g., oxygen vacancies); yet asymptotically a truly superconducting 
low-temperature "vortex glass" phase has been predicted [3] . 

Therefore, in addition to point defects, the influence of extended or corre- 
lated disorder, promising stronger pinning effects, on vortex transport proper- 
ties has been considered. Experimentally, linear damage tracks have been pro- 
duced in oxide superconductors by irradiation with high-energy ions. These 
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columnar defects serve as strong pinning centers for the flux lines, thus sig- 
nificantly increasing the critical current and shifting the irreversibility line 
upwards (for references, see Ref. [1], Sec. IX B). 

In a long-wavelength description in the spirit of the London approxi- 
mation, one may consider the following model free energy for N flux lines, 
defined by their trajectories r^z) as they traverse the sample, interacting 
with each other and with Ad columnar defects which are aligned along the 
magnetic-field direction z [4], [5]: 



Here rij(z) — \ri(z) — rj(z)\, and V(r) = 2e Ao(r / A) denotes the repulsive 
vortex interaction; the modified Bessel function K (x) cx — lnx for x — ► 0, 
and cx x~ x / 2 e~ x for x — > oo. Thus the London penetration depth A defines 
the interaction range and the energy scale eo = {4>o/^ A) 2 . For low fields 
B < <po/^ 2 , i- e - : A iS a = (4/3) 1 / 4 (^o/5) 1 ^ 2 , the vortex line tension is ei w e 
(a is the triangular lattice constant in the corresponding pure system). Fi- 
nally, the pinning energy is a sum of Ad z-independent potential wells Vb (r) 
with average spacing d centered on the positions {Rfe}, whose diameters are 
typically bo « 100 A, with a variation of Sbk/b n w 15%, caused by the ion- 
beam dispersion. This induces some distribution P of the pinning energies 
Uk, which may be determined from the formula Uk ~ (eo/2) ln[l + (6fc/\/2£) 2 ]i 
where £ denotes the coherence length [5]. The thermodynamic properties of 
the model (1) are to be obtained by evaluation of the grand-canonical parti- 
tion function 



with the chemical potential fx = (H Cl — H)<Po/4tt, which changes sign at the 
lower critical field H Cl = 4>o ln(A/^)/47rA 2 . 

As for pure systems [2] , one can map this statistical problem of directed 
flux lines interacting with columnar pins, onto the quantum mechanics of 
bosons in two dimensions subject to static point disorder via the identification 
of the vortex trajectories r^z) with the particle world lines in imaginary 
time [4], [5], [6]. In Table 1, the corresponding quantities for bosons and 
magnetic flux lines are listed (notice that the thermodynamic limit L — > 
oo for the directed polymers corresponds to a zero- "temperature" quantum 
problem), and the ensuing phase diagram is sketched schematically in Fig. 1. 
At high temperatures, one thus finds an entangled liquid of unbound flux lines 
(corresponding to a boson superfluid), separated by a sharp second-order 
transition (sec Ref. [7] and references therein) from a low-temperature phase 
of localized vortices. This Bose glass phase is characterized by an infinite 
tilt modulus C44, and turns out to be stable over a certain range of tipping 
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angles of B away from the z direction [5] , [8] . At least for A <C d, the theory 
also suggests a low-temperature Mott insula tor phase, when B = B^, i.e., the 
vortex density exactly matches that of the columnar pins. The Bose glass and 
the Mott insulator are distinct thermodynamic phases, for the latter should 
be characterized by infinite tilt and compressional moduli [5] . 



Table 1. Boson analogy applied to vortex transport. 

Bosons m h /3h V(r) n p p B E j 

Vortices ei T L 2e K (r/\) B/fo (H — H Cl )4>o/^ p 2 c^ z x J/c £ 



Fig. 1. Schematic phase diagram for vortices interacting with columnar defects 
aligned along the magnetic-field direction. 



The formal analogy with a boson superfluid can be further exploited to 
investigate the effect of disorder on density and current correlation functions 
in the flux liquid phase within the Bogoliubov approximation [9]. Thus one 
may compute the corresponding transport coefficients, e.g., the boson su- 
perfluid density p s whose flux line analog is the inverse tilt modulus of the 
vortex array (see Table f ). For parallel columnar defects of strength Uq one 
eventually finds [10]: 

, withT B G(B)«T*[ ( /»o/(47rA) 2 B] 1 / 4 (3) 

and k^T* = b(iiUo) 1 ^ 2 , in accord with estimates obtained in the Bose glass 
itself [5] . The critical behavior at the boson localization transition is not yet 



c^(T) « (ner)" 1 1 - (T BG /T) 4 
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fully understood, although quantum Monte Carlo simulations [11] seem to 
support the scaling theory [7] . Direct simulations for the critical dynamics of 
flux lines near the Bose glass transition have also been performed [12]. 

2 Structural Properties in the Bose Glass Phase 

Recently, the positions of flux lines and parallel columnar pins, aligned along 
the magnetic-field direction, were simultaneously measured in a BSCCO sam- 
ple in the Bose glass phase at low magnetic fields using a combined chemical 
etching and magnetic decoration technique [13]. Figure 2 depicts the thus ob- 
tained positions of TVd columnar defects and N flux lines in a two-dimensional 
cross section perpendicular to B in a case where / = N/N^ = B/B^ rts 0.24 
and A w 0.45 a . While the columnar defects are to a good approximation 
randomly distributed in space, the flux lines clearly form a highly correlated 
configuration; correspondingly, the two-dimensional vortex structure func- 
tion S(q) displays a distinct peak at wave vector q w 2n/a a , resembling an 
amorphous material (see Fig. 3). This correlated structure must obviously be 
the result of the mutual repulsive interactions Vint between the vortices which 
dominate over any statistical fluctuations in the pinning potentials 5Vr>; on 
the other hand, all the flux lines are attached to a defect, and therefore 
Vu >V int [14]. 



Fig. 2. Positions of empty columnar defects (open circles), and pins occupied by 
flux lines (filled circles), as obtained from a decoration experiment (/ « 0.24). 



In order to model this experimental situation, we can fortunately simplify 
the free energy (1) substantially: (i) As long as we limit ourselves to "low" 
temperatures (an estimate for BSCCO actually yields T T\ = 0.9 T c ) 
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we can neglect thermal wandering of the flux lines, which is furthermore 
significantly suppressed by the confining line defects, we may drop the tilt- 
energy oc ei and consider an effectively "classical" (ksT ~ h in the boson 
analogy) two-dimensional system, (ii) For magnetic fields smaller than the 
crossover field B*{T) in Fig. f, the vortices are sufficiently dilute to ensure 
that their repulsive interactions have no influence on the localization length, 
and in fact the "boson" statistics becomes irrelevant (for BSCCO, B*(T\) w 
0.7 B<f,). Having thus disposed of the entropic corrections, we may consider the 
following effective Hamiltonian and its grand-canonical counterpart [14], [15], 



where i = 1, . . . , N-£, denote the defect sites, randomly distributed in the x- 
y plane, = 0,1 is the corresponding site occupation number, and the ij 
are random-site energies (originating in the varying pin diameters), whose 
distribution of mean zero and width w w 0.1/eo we for simplicity choose to 
be flat: P(U) = (w - /2w. 

This classical model is still far from trivial, however, due to the inter- 
play of long-range interactions and disorder. Fortunately models of the form 
(4), albeit with a Coulomb potential replacing the screened logarithmic in- 
teraction V(r), have been extensively studied in the context of disordered 
semiconductors, which allows us to adapt qualitative arguments and numeri- 
cal simulation procedures from the literature (see, e.g., Refs. [16], [17]). Thus, 
in order to infer the (approximate) ground-state properties in the Bose glass 
phase as functions of the filling / and interaction range X/d we employ a 
zero-temperature Monte Carlo algorithm that essentially minimizes the to- 
tal energy (4) with respect to all possible one-vortex transfers; i.e., the ground 
state stability is checked by demanding that all single-particle hops from oc- 
cupied sites i to empty sites j require an energy Zl,^ = ej — ei — V(rij) > 0, 
where = U + n jV{ r ij ) are the interacting site energies. For A <J d 

and not too large / 0.6 the resulting distribution of flux lines among the 
columnar pins closely resembles the experimental data in Fig. 2. More quan- 
titatively, we can investigate the amorphous peak at qcio ~ 2tt in the vortex 
structure factor 5(q) = J2i j e iq '( r<_r ^/iV, averaged over all directions in re- 
ciprocal space. Whereas the dependence on the interaction range is rather 
weak as long as A J2 d, the peak in Fig. 3 vanishes when / *Z 0.4; this reflects 
that at higher fillings the system has to increasingly accomodate with the 
underlying randomness, and therefore the correlations induced by the inter- 
actions disappear. Furthermore, a triangulation of the vortex positions shows 
that typically only about 50 % of the sites are sixfold coordinated, and thus 
the structure formed by the flux lines rather resembles an amorphous state 
(justifying the name "Bose glass") rather than a weakly disordered triangular 
lattice [15]. 
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Fig. 3. Dependence of the vortex structure factor peak on the filling / for X/d = 2; 
solid: / = 0.1, dot-dashed: / = 0.2, dashed: / = 0.4, and long-dashed: / = 0.8. 

3 The Distribution of Pinning Energies 

The above Monte Carlo algorithm can also be employed to determine the 
distribution of pinning energies g(e), with the vortex interactions taken into 
account. As in the boson picture this is nothing but the single-particle density 
of states, g(e) is readily obtained by sampling the site energies for a number 
of disorder realizations and then performing a quenched average. Following 
the mean-field arguments of Efros and Shklovskii (1984) for a system with 
long-range repulsive interactions V(r) cx r~ CT , with < a < D (here, the 
space dimension is D = 2), one expects that the density of states vanishes at 
the chemical potential according to a power law, 

5 (e) cx |e- M | Srff , (5) 

with a gap exponent smf = (D/<r) — 1. This means that as a consequence 
of the correlations induced by the long-range interactions, the probability 
of finding a low-energy empty site (e iJ fi) for introducing an additional 
particle is greatly reduced, as is the density of states of high-energy filled 
sites (e fi) for the corresponding process of removing a particle from the 
system. While the qualitative prediction (5) remains correct, the actual gap 
exponent is influenced by subtle fluctuation effects [17] and therefore differs 
from the mean-field predictions (for the currently most accurate results for 
Coulomb interactions, a = 1, see Ref. [18]). 

Figure 4 depicts two characteristic examples for the distribution of pin- 
ning energies <?(e), as obtained from our simulations [15], and normalized 
according to J g(e)de = l/d 2 . In the case of an infinite-range logarithmic po- 
tential, the Ewald summation method was used, and before taking the limit 
A — > oo, a term e- — ej = —2-Kf(X/d) was subtracted from the site energies. 
The filled and empty states are separated by a wide soft gap centered at the 
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Fig. 4. Normalized distribution of pinning energies G(E) — 2e d 2 g(e) as function of 
the single-particle energies E = e/2eo, with (a) A — > oo, and / = 0.2, (b) \/d = 2, 
and / = 0.4. The location of p is marked by the arrow. 

chemical potential, with an effective gap exponent s e s w 3 for low filling. 
Remarkably, this characteristic double-peak structure and the marked de- 
pletion of states persists even for finite-range repulsion with A > d (Fig. 4b), 
and although g(p) does not strictly vanish any more, the distribution of pin- 
ning energies for e w p can be fitted with a power law (5); e.g., s e ff ~ 1-2 
for X/d = 1 and / = 0.1. Upon increasing /, the gap closes quickly due to 
the stronger influence of the underlying randomness (compare Fig. 3), but 
for / ^ 0.2 the correlation effects only disappear when A < d. 

4 Variable— Range Hopping Transport of Flux Lines 

An in-plane current JIB induces a Lorentz force = x J/c, acting on 
the vortices. Accordingly, a term 5Fn = — dz ^ i f l • r, (z) has to be added 
to the free energy (1). In the spirit of the thermally assisted flux-flow (TAFF) 
model of vortex transport, we write the superconducting resistivity (i.e.: the 
conductivity in the boson representation) as p = £/ 'J w pocxp [—Ub(J)/T], 
where po is a characteristic flux-flow resistivity, and Ub represents an effective 
barrier height. We shall see that in the Bose glass phase Ub{J) actually 
diverges as J — > 0, 

£ w p Jcxp{-E Q /kBT)(J /JT^} . (6) 

Driven by an intermediate external current J\ < J < J c , a flux line will 
start to leave its columnar pin by detaching a segment of length z into the 
defect-free region, thereby forming a half-loop of transverse size r (Fig. 5a). 
The free energy corresponding to this saddle-point configuration is readily 
estimated to be 8F\{r, z) ss iir 2 / z + Uqz — j\jz, and optimizing successively 



8 Uwe C. Tauber 



Fig. 5. Saddle-point configurations for vortex transport in the Bose glass phase, 
(a) Vortex half-loop, (b) Double-superkink in the variable-range hopping regime. 



with respect to z and r one finds z* w r* y^ei/Ua s=s c^/TJJq/^qJ . Hence 
the current -voltage characteristics assumes the form (6), with Eq = Ek = 
dyeiUo, Jo = Ji ~ cUo/(f>od, and the effective transport exponent p c s = 1 in 
the vortex half-loop regime [5] . For J — > J\ , the flux line will extend a double 
kink of width wk — d^ii/Uo and energy 2Ek to an adjacent columnar defect, 
and therefore for J < J\ neighboring pins must be taken into account. 

Thus, at low currents one enters another regime, where the flux line emits 
a pair of kinks to a possibly distant defect which happens to provide a fa- 
vorable pinning energy (Fig. 5b). This is the vortex analog of variable-range 
charge transport in disordered semiconductors (see Ref. [16] and references 
therein). The cost in free energy for such a configuration of transverse size R 
and extension Z along the magnetic-field direction will consist of the double- 
kink energy 2Ek(R) — 2E-^{d)R/d stemming from the elastic term, and the 
difference in pinning energies of the highest-energy occupied site, et w fj,, and 
the empty site at distance R with ej = /i + A(R), where the distribution of 

pinning energies enters through the condition j^ +A ^ g(e)de = R~ D . Min- 
imizing the saddle-point free energy <5Fsk ~ 2E\^R/d+ ZA(R) — JlRZ in 
the regime where g(e) has the form (5), one arrives at (6) with E = 2Ek 
and the effective transport exponent [15] 

PeS = (Scff + 1)/(D + Scff + 1) • (7) 

For short-range interactions, g(p) > 0, i.e., s ff = 0, and the Mott 
variable-range hopping exponent pu = l/(D + 1) is recovered [5]. But as 
the interaction range A becomes larger than the average vortex separation 
ao = d/y/f, the emerging correlations reduce the phase space for single- 
vortex hopping drastically, and p c g < pm', in addition, as the normalized 
distribution of pinning energies is broadened by the interactions, also the 
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Fig. 6. Log-log plots (base 10) of the activation energy Ub{J) oc R*(J)/d for vari- 
able-range hopping vs the reduced current j = J(f>od/2eoc for / = 0.1. Diamonds: 
A — > 0, circles: X/d = 1, squares: X/d = 2, triangles: X/d = 5. 

current scale Jo in (6) is increased. The result of the minimization proce- 
dure using the numerical data for g(e) in the case / = 0.1 is shown in Fig. 6 
as a function of the interaction range [15]. Assuming that flux line trans- 
port in the low-current regime is still dominated by single-vortex hops, this 
implies that long-range repulsive interactions reduce the resistivity by sev- 
eral orders of magnitudes with respect to the non-interacting situation. E.g., 
for log j w 1.5 the exponent in (6) is about ten times smaller for X/d = 5 
as compared to A — > (see Fig. 6), and hence p is reduced by a factor of 
w 10~ 5 . In accord with the relation (7), the effective transport exponent as- 
sumes values upto p c g rj 0.7 for A — > oo and / = 0.1, while, e.g., for X/d = 1 
and / = 0.1 we find p c ff ~ 0.55. We remark that the former value for for 
the variable-range hopping exponent of logarithmically interacting particles 
in two dimensions is consistent with an analysis based on the very different 
assertion that actually collective many-particle hops yield the lowest energy 
barriers [19]. Finally, a recent experiment on a BSCCO sample with / rj 0.15 
and X/d « 1.6 reported a measured p c s 0.57 in the variable-range hop- 
ping regime across about half a decade in the current -voltage characteristics 
[20] ; this apparently agrees well with our results, although due caution needs 
to be applied for this comparison as in this experiment the magnetic field 
distribution was in fact very inhomogeneous. Yet further experiments should 
clearly be capable to test the above predictions quantitatively. 

In summary, experiments have established that columnar defects provide 
very effective pinning centers for flux lines in high-T c superconductors, and 
demonstrated that in certain cases the repulsive vortex interactions lead to 
remarkable spatial correlations. As a consequence, a soft "Coulomb gap" 
emerges in the distribution of pinning energies, and vortex transport in the 
variable-range hopping regime should be drastically reduced. 
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